A highly-collimated SiO jet in the HH212 protostellar outflow 
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ABSTRACT 

We mapped the inner 40 of the HH212 Class outflow in SiO(2-l), SiO(5-4) and continuum 
using the Plateau de Bure interferometer in its extended configurations. The unprecedented 
angular resolution (down to 0'.'34) allows accurate comparison with a new, deep H2 image obtained 
at the VLT. The SiO emission is confined to a highly-collimated bipolar jet (width ~ 0'.'35) along 
the outflow axis. The jet can be traced down to within 500 AU of the protostar, in a region that 
is heavily obscured in H2 images. Where both species are detected, SiO shows the same overall 
kinematics and structure as H2, indicating that both molecules are tracing the same material. We 
find that the high-velocity SiO gas is not tracing a wide-angle wind but is already confined to a 
flow inside a narrow cone of half-opening angle <6° at < 500 AU from the protostar. Transverse 
cuts reveal no velocity gradient compatible with jet rotation above 1 km s" 1 , in contrast to 
previous claims based on H2 spectra. 

Subject headings: Stars: formation - Radio lines: ISM - ISM: jets and outflows - ISM: molecules - ISM: 
individual objects: HH212 



1 



1. Introduction 

The launching of jets from young stars is one 
of the most enigmatic and intriguing phenom- 
ena in astrophysics. One fundamental problem to 
which jets are believed to provide a solution is the 
removal of angular momentum from the central 
accreting system, which is required to permit low 
angular momentum material to f all onto the star 



Rece nt optical observations (e.g. iBacciotti et al 



|2002( ) suggest that atomic jets from evolved T 
Tauri stars may indeed transport angular momen- 
tum away from the inner regions of the accre- 
tion disk. Is the same mechanism at the origin 
of the molecular jet counterparts observed when 
the protostar is still deeply embedded in its natal 
high-density core? To address this question, we 
have investigate HH212 in Orio n (d = 460 pc), a 



highl y symmetric bipolar H2 jet (jZinnecker et al 



1998), associated with a collimated CO outflow 



( Lee et al. 20061) , and driven by a low-luminosity 
Class protostar, IRAS05413-0104. HH212 pro- 
vides an optimal situation in which to study 
jet kinematics, si nce it lies close to t he plane 
of the sky (4°; Claussen et al. 19981) . shows 



hints of rotation in one H2 knot (jDavis et al. 
2000), and is surrounded by a compact ammo- 



nia core rotating about an axis aligned with the 
jet (jWiseman et al. 200lh . As a tracer, we used 
SiO, which is associated with shocks and gener- 
ally suffers minimal contaminat ion from infalling 



envelopes or swept-up ca vities ([Guilloteau et al. 
Il992t iHirano et al. Il200fih . 



2. Observations 

The SiO observations of HH212 were obtained 
in January-March 2006 with the IRAM Plateau 
de Bure Interferometer (PdBI) in France. Two 
6-element configurations were used: the new ex- 
tended A configuration, which includes baselines 
from 136 m up to 760 m, and the new B con- 
figuration, with baselines from 88 m to 452 m. 
The dual-channel receivers and the correlator were 
tuned to simultaneously observe the SiO J = 2-1 
line at 86.847 GHz (resolution: ~ 0.13 km s _1 ; 
bandwidth: ~ 70 km s _1 ) and the J = 5-4 line 
at 217.105 GHz (0.11 and 55 km s -1 ). Continuum 
emission at both frequencies was also measured. 
Amplitude and phase were calibrated by observing 
0528+134 and 0605-085. The weather conditions 



were excellent (phase noise rms of 30° at 86 GHz 
on the 760 m baselines), leading to a radio seeing 
of 0'.'35. The flux density scale was derived by ob- 
serving 3C 84 and 3C 279, with an uncertainty of ~ 
25%. Images were produced using natural weight- 
ing, and restored with clean beams of l'/89 x 0'.'94 
(PA=22°) at 3.5 mm, and of 0'.'78 xO'/34 (PA=21°) 
at 1.4 mm, using the GILDAS software. The high 
degree of elongation of the beam is due to the zero 
declination of the source. However, the spatial 
resolution is highest perpendicular to the jet di- 
rection, making the present observations ideal for 
the study of the jet collimation and kinematics. 

A new, deeper, and higher spatial resolution 
NIR image of HH212 was obtained using ISAAC 
on the ESO Very Large Telescope (VLT), Paranal, 
Chile on November 12, 2005 under good condi- 
tions. Although the knots have proper motion s 



of 100-200 km s" 1 (jMcCaughrean et al. 1 12003) . 



the small epoch difference to the SiO data en- 
sures that the knots have moved only ~ 0'.'025, 
thus facilitating an comparison between the two 
data sets. A series of images were taken through 
a 1% wide filter centred on the v=l-0 S(l) line 
of H2 at 2.122/im. with a total integration time 
of 32 minutes per pixel after mosaicing. The im- 
age scale was 0'.'147 pixel -1 , and the spatial res- 
olution is 0'.'35 FWHM, similar to our SiO (5-4) 
map. Accurate astrometry (0'.'14 rms) was de- 
rived via 21 stars in common with the 2MASS 
Point Source Catalog. The images have not been 
continuum-subtracted and thus stellar sources are 
visible. However, the HH212 jet itself consists al- 
most exclusively of shocked H2 line emission and 
thus there is no confusion along the jet axis. 

3. Results and Discussion 

Emission maps of SiO J = 2-1 and 5-4 are 
shown in Fig. la-b, superimposed onto the nearly 
contemporaneous H2 image. The SiO emission is 
confined to a highly collimated bipolar jet, located 
along the main outflow axis. The red- and blue- 
shifted lobes consist of two symmetric pairs of bul- 
lets: an outer pair (Blue 1 and Red 1) that follows 
the H2 intensity distribution beyond 10" from the 
driving source, and a new inner pair within ±2" of 
the source (Blue 2 and Red 2), with no H2 coun- 
terpart. Figure lc shows the inner 5" of the 1.4 
mm continuum map, with the SiO (5-4) emission 
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superimposed. 

3.1. Continuum sources 

The 1.4 mm continuum map (Fig. lc) shows a 
bright source (hereafter called MM1) at position 
a J20 oo= 05 h 43 m 51'41, 5 J20 oo= -01° 02' 53'/160, 
in excellent agree ment with the VLA pos i tion a t 
3.5 cm reported by IGalvan-Madrid et al. I (|2004h . 
Comparison with the SiO(5-4) map clearly indi- 
cates that MM1 is the driving source of the molec- 
ular jet. In addition, we have a tentative detec- 
tion at S/N ~ 6 (~ 4 mJy) of a secondary peak 
(hereafter MM2) at Aa=+l'/6, A5=-CX'5, which 
could trace an embedded companion. Only MM1 
is bright enough to be detected at 3.5 mm, with 
a flux of 6 mJy. Gaussian fitting in the UV plane 
at 1.4 mm shows that MM1 is associated with an 
unresolved peak (< 0'.'26) with a faint extension 
towards SE. The integrated flux at 1.4 mm is 33 
mJ y, significan t ly low er than the 110 mJy reported 
by iLee et al. 



~ 0'.'35, and are thus much more collimated than 
the hourglass-shape d bipolar cavity se en in 13 CO 



on similar scales by Lee et al. ( 20061 ). Position- 



Velocity (PV) diagrams along the jet axis (Fig. 2) 
further show that the SiO emission extends to high 
velocities of ±1 km s -1 from the amb ient velocity 



(+1.6 km s 1 ; I Wiseman et al. I | 200lh. compare d 



to only ± 1.7 km s" 1 for 13 CO (|Lee et al. 1120061 ) . 
Hence the SiO lobes include a narrower and faster 
jet-like component distinct from the swept-up cav- 
ity. The maximum radial velocity is similar to the 
typical centroid velocity of H2 knots further out 
(see Fig. 2), arguing that the high- velocity SiO is 
probably tracing the base of the large-scale (~ 0.5 
pc) molecular jet. Co rrection for inclination (4° to 
the plane of the sky; IClaussen et al. 19981 ) gives 
a full jet speed ~ 140 km s " 1 , comparable to the 
prop er motion of H2 knots (jMcCaughrean et al 



2002). Inner SiO knots were also recently im- 



aged in HH211 down to < 2" from the p rotostar 



,iy luwui Ulan Luc uumjj lepmueu r~ 11 1 i = m -r 

d2006h in a 2'.'5 beam, indicating (|Hirano et al. || 2006 fc | Gueth et al. | |2006j ). Hence 



that most of the envelope emission has been re- 
solved out by the PdBI. Contamination by free- 
free emission from the jet is negligible: if we ex- 
trapolate the flux measured b y the VLA at 3.5 cm 



(jGalvan-Madrid et al. I l2004h with a typical v 



,0.6 



slope, the expected contribution is only 0.16 mJy 
at 3.5 mm. Hence we are dominated by dust con- 
tinuum. The 1.4 mm flux corresponds to a mean 
surface brightness of 3.4 K in our beam. Since 
envelopes and disks are not elongated along the 
jet axis, the source size is less than 0'.'26 x 0'.'26 
and thus the true brightness temperature is quite 
high: > 13 K. The spectral slope between 1.4 mm 
and 3.5 mm, derived from UV data, is ~ 2. Both 
facts seem to indicate that dust emission is close 
to optically thick at 1.4 mm. This would suggest 
that the unresolved mm peak is not tracing the in- 
ner parts of the envelope, but a circumstellar disk 
viewed close to edge-on. The disk diameter would 
be < 0'.'26 = 117 AU. 

3.2. The inner SiO bipolar jet 

Our high resolution maps reveal a small-scale 
pair of SiO emission knots emerging from MM1 
and extending out to l"-2" = 500-1000 AU 
(Fig. lc). No H2 counterpart is visible, as the 
whole structure lies entirely in the high extinction 
region around the source. The beam-deconvolved 
transversal size of the two lobes (Blue 2/Red 2) is 



SiO appears to be a powerful tracer of the jet base 
in Class sources. 

Hirano et al. ( 20061 ) found broad SiO linewidths 
in the protostar in HH211, which they attributed 
to a wide angle wind. Such an interpretation is not 
supported in HH212, however: SiO profiles show 
little blue/red overlap in each lobe (see Fig. 2), de- 
spite the almost edge-on inclination. Assuming a 
conical wind of constant speed, a half-opening an- 
gle < 6° is inferred, ruling out a wide-angle wind. 
The fact that SiO emission stops 'precisely near 
zero velocity, with an overlap of at most 2 km s , 
would actually suggest that the wind half-opening 
angle coincides to within 2° with the angle from 
the plane of the sk}0. Since such a coincidence 
is statistically unlikely, it appears more probable 
that high-velocity SiO is confined to less than 4° 
from the jet axis (thus producing no emission at 
zero velocity) and that SiO emission near systemic 
velocity is tracing intrinsically slower, possibly less 
collimated material. This slow material may trace 
unresolved bowshock wings. Indeed, water maser 
spots observed via VLBI reveal curved bowshocks 



1 The wind half-opening angle, ma x, is related to the angle 
to the plane of the sky, a, by tan6> max = (R — 1)/(R + 
1) tan a where R is the ratio (in algebraic value) of max- 
imum to minimum radial velocities in a given lobe. In 
HH212, ti max > 10 km s _1 , -2 < v min < km s" 1 , and 
a=4° yield R < -5 and max = 4°-6°. 
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at ± O'.'l from the protostar, covering a range of 
radial v elocities (black arrows in the right panel of 



Fig. 2; Claussen et al. 1998f ). Given the restric- 



tive excitation and coherence conditions for maser 
amplification, the observed velocity range is less 
than the full range in the bow. Hence, internal 
bowshocks could contribute significantly to the 
broad SiO line widths. Intrinsic gradients in jet 
speed (due e.g. to time variability, or to a range 
of launch radii) may also be present. 

3.3. The outer SiO bullets 

At distances of 10" and beyond, the spatial cor- 
respondence between SiO and H2 is very good. 
The two symmetric elongated SiO bullets (Blue 
1/Red 1) peak towards H 2 knots SK2-SK4 and 
NK2-NK4 (note that the SiO (5-4) image is af- 
fected by primary beam attenuation at large dis- 
tances from the centre (FWHM = 22"). SiO emis- 
sion is also seen between and beyond these knots, 
with a morphology similar to that in H 2 (Fig. 
1). The correspondence between SiO and H 2 in 
the outer bullets is also excellent in radial veloc- 
ity, as shown in the PV diagram in Fig. 2 (left 
panel). The SiO (2-1) profiles peak within 2 km 
s _1 of the H2 centroid. Hence SiO appears to trace 
the same jet a s H 2 , confirming the conclusion of 
Takami et al~l (|2006f) based on SiO (8-7) profiles 
at much lower angular resolution (22"). The SiO 
line profiles are broader than towards the inner 
knots, but they show the same extension down to 
zero velocity and small degree of red/blue over- 
lap, suggestive of internal bowshocks. Line profile 
asymmetries between the two lobes indicate slight 
differences in shock structure despite the highly 
symmetric knot spacings. 

In contrast to the excellent agreement between 
SiO and H 2 beyond 10", we find a striking lack of 
SiO emission towards the brightest H 2 knots (SKI, 
NK1). A similar effect was observed towards the 
HH211 jet, where the brightest, b ow-shaped H 2 



knots also lack an SiO counterpart (jHirano et al 
200G). One possible explanation would be that 



SKI and NK1 trace more powerful shocks where 
SiO does not form/survive or where the SiO ex- 
citation is extremely high and consequently the 
low-J transitions are very weak. The detection of 



[Fe II ] emission in SKI and NK1 (Zi nnecker et al. 



19981: ICaratti o Garatti et al~l 120061 ). their clear 
bow-shock geometry (see Fig. 1), and their 3-4 



times broader linewidth compared to other H 2 
knots in HH212 (Fig. 2) support this hypothesis. 
Estimates of SiO excitation and abundances will 
be presented in a separate forthcoming paper. 

3.4. Search for jet rotation 



Davis et al. 



(|2000h reported a centroid veloc- 
ity shift of +2.3 km s _1 from NW to SE across 
knot SKI which they attributed to jet rotation, 
based on its agreement with th e rotation sense 



of the ammonia protostellar core (jWiseman et al. 
120011 ). Although SKI is not detected in SiO (see 
3.3), we can test this interpretation in the other 
knots, as rotation speeds should remain roughly 
constant onc e the jet has achieved cylindrical col- 



limation (see lPesenti et al. 1120041) 



Figure 3a presents PV diagrams across the jet 
in SiO (5-4) for each emitting region. No shift of 
the intensity peak is seen a cross Blue 1 on th e 
scales ± 0"5 investigated by iDavis et al. ( 2000h . 
On the other hand, a clear gradient is observed 
across the Red 1 bullet (Fig. 3, bottom panel), 
with a shift of -6 km s _1 over ~ 1" from NW to 
SE. However, the gradient goes in opposite sense 
to the rotation pattern of the NH3 core, ruling out 
that we are observing rotation of a centrifugally- 
driven MHD jet. The fact that the centroid of 
Red 1 is slightly displaced from the jet axis, i.e. 
that outer SiO bullets are slightly misaligned with 
respect to the inner SiO jet, suggests that this 
effect is due instead to a kinematical asymmetry 
in the shocked material, possibly related to the 
H 2 jet wiggling apparent in Fig. 3b. Since SKI 
traces a much wider bowshock than the SK2-SK4 
knots in Red 1 (see Fig. 1), its transverse gradient 
measured in H 2 could easily be affected by similar 
asymmetries, and should be interpreted with great 
caution. The innermost SiO knots (Blue 2/Red 
2) do not show signs of transverse velocity shifts 
above 1 km s _1 , consistent with the high- velocity 
SiO jet being practically unresolved transversally. 
Higher angular resolution would be needed to con- 
strain the level of jet rotation. 
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Fig. 1. — a: Contour map of integrated emission in SiO(2-l) superimposed on the new H2 image. The SiO 
emission regions are labelled in white. The filled triangle at (0",0") denotes the driving source MM1 (see 
Sect. 3.1). Contour levels range from 3a (0.22 Jy beam -1 km s" 1 ) to 39<r by steps of 6a. The filled ellipse 
in the lower left corner shows the synthesised PdBI beam (HPBW): l'/89 x 0'.'94. b: same as panel (a) for 
the SiO (5-4) emission. The dotted circle indicates the primary beam of the PdBI antennae. H 2 knots are 
labelled in black. Contour levels range from 3a (0.42 Jy beam" 1 km s" 1 ) to 39cr by steps of 3a. The beam 
is 0"78 x 0'.'34. c: Inner 5" of the SiO(5-4) map, superimposed onto the continuum image at 1.4 mm (gray 
scale and white contours). Continuum contour levels range from 3a (1.9 mJy beam" 1 ) to 36a by steps of 6a. 
Labels indicate the HH212 driving source (MM1) and a tentative weaker source (MM2).[5ee the electronic 
edition of the Journal for a colour version of this figure.] 
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HH212 SiO(2-1) outflow main axis Si0(5-4) 
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Fig. 2. — Left: Position-velocity cut of SiO J=2-l along the jet axis (PA = 22°). Contour levels range from 
1.4 K to 27.5 K by 2.9 K (8a). Dashed lines mark the position of MM1 and the cloud VLsr (+1-6 km s _1 ; 
Wiseman et al. 2001). Filled triangles with error bars show the velocity centroid and intrinsic FWHM of 
H2 knots (Takami et al. 2006). Right: Zoomed- in PV cut of the inner jet in SiO J=5-4. Contour levels 
range from 3 K to 21 K by 3 K (4cr). Black arrows show the velocity range of H 2 masers ±0.1" from the 
protostar (Claussen et al. 1998). 
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SiO(5-4) H 2 1-0 S(1 ) 




Angular Offset (arcsec) 

Fig. 3. — a: PV cuts in SiO(5-4) perpendicular to the jet axis, averaged over each of the four SiO bullets. 
Contour levels start at 3cr with steps of 6a. They range from 1 K to 9.3 K (Blue 1), 2.7 to 14.5 K (Blue 
2/Red 2), and 1.7 to 11.6 K (Red 1). Dashed lines mark the jet axis at PA=22° and the VLsr °f the cloud 
(+1.6 km s _1 ; Wiseman et al. 2001). A solid line shows the velocity gradient across Red 1. b: Zoom-in of 
the H 2 image highlighting the jet wiggling. Vertical black arrows indicate the declination range of the Blue 
2/Red 2 inner jet. 
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